NONLINEAR SPIN DYNAMICS IN NUCLEAR MAGNETS 



V.I.Yukalov 

Department of Mathematics and Statistics 
Queen's University, Kingston, Ontario K7L 3N6, Canada 

and 

Bogolubov Laboratory of Theoretical Physics 
Joint Institute for Nuclear Research, Dubna 141980, Russia 



PACS: 76.20+q; 76.60Es 



1 



Abstract 

A method is developed for solving nonlinear systems of differential, or inte- 
grodifferential, equations with stochastic fields. The method makes it possible to 
give an accurate solution for an interesting physical problem: What are the pe- 
culiarities of nonlinear spin dynamics in nonequilibrium nuclear magnets coupled 
with a resonator? Evolution equations for nuclear spins are derived basing on a 
Hamiltonian with dipole interactions. The ensemble of spins is coupled with a 
resonator electric circuit. Seven types of main relaxation regimes are found: free 
induction, collective induction, free relaxation, collective relaxation, weak super- 
radiance, pure superradiance, and triggered superradiance. The initial motion of 
spins can be originated by two reasons, either by an imposed initial coherence or 
by local spin fluctuations due to nonsecular dipole interactions. The relaxation 
regimes caused by the second reason cannot be described by the Bloch equations. 
Numerical estimates show good agreement with experiment. 
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I. Introduction 

Spin dynamics in polarized nonequilibrium systems is usually described by using 
the Bloch equations for the components of uniform magnetization. The derivation of 
the Bloch equations from the evolution equations for a spin model can be found, for 
example, in ter Haar [1]. The solution of the Bloch equations for the case of small 
deviations from a stationary state is straightforward and well known in the theory of 
magnetic resonance [2]. The situation becomes more complicated when the spin system 
is coupled with a resonator. Then there appears an essential nonlinearity due to the 
action of resonator feedback field. The system of the coupled Bloch and resonator-field 
equations is typical of the theory of maser amplifiers and generators [3]. 

The nonlinear system of the Bloch and resonator-field equations can be slightly 
simplified by invoking the slowly-varying amplitude approximation [3]. However, this 
does not help much, since the resulting equations are, as before, nonlinear. To achieve 
further simplification, one resorts to the adiabatic approximation which leads to the 
proportionality of the feedback field to transverse magnetization, that is, to the static 
coupling [4-7]. The adiabatic approximation, as is known [8], works well only at the final 
stage of relaxation processes when different variables adiabatically follow each other, 
but it cannot correctly describe intermediate stages where transient phenomena occur. 
The incorrectness of the static-coupling approximation is physically evident, as only 
moving spins, but not immovable, are able to induce a field in resonator. More accurate 
is the dynamic-coupling approximation [9] in which the feedback field is proportional 
to the time derivative of transverse magnetization. But both these, static-as well as 
dynamic-coupling, approximations do not take into account retardation effects that 
may be important for transient phenomena. 

Moreover, the Bloch equations themselves may be inappropriate for explaining some 
kinds of relaxation processes. This concerns, for example, the interpretation of the 
recent series of experiments [10-15] observing nuclear spin superradiance. In these 
experiments a nonequilibrium system of polarized nuclear spins is placed inside a coil of 
a resonance electric circuit. The initial polarization is directed opposite to an external 
magnetic field. If this polarization is sufficiently high and the coupling with a resonator 
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is enough strong, then the power of current, as a function of time, after some delay, 
displays a sharp burst with a damping time much shorter than the dephasing time 
T 2 . This time behaviour of the current power is analogous to that of the radiation 
intensity of atoms or molecules in the case of optical superradiance. Because of this 
analogy, the corresponding coherent phenomenon in spin systems has also been called 
superradiance, or more concretely, spin superradiance. Friedberg and Hartmann [16] 
pointed out that the whole process of interaction of a spin system and a resonance coil, 
in fact, involves no radiation into free space but merely nonradiative transfer of energy 
from the sample to the coil, where the energy is dissipated ohmically. Nevertheless, 
the term spin superradiance has become commonly used. The excuse for this is not 
solely the formal analogy of temporal behaviour of current power, for spin systems, 
and of radiation intensity, for atomic and molecular systems, but also a deep physical 
similarity: The spin superradiance, as well as optical superradiance, is a collective 
process of coherent self-organization. Although the self-organized coherence of spin 
motion develops not because of a common radiation field, as in atomic and molecular 
systems, but owing to a resonator feedback field. In addition, coherent motion of 
spins inevitably produces coherent magnetodipole emission with properties completely 
analogous to superradiance of optical systems, though the magnetodipole radiation 
intensity is too weak to be measured as easy as the power of current [17]. 

In the same way as for optical systems [18], one has to distinguish the pure from 
triggered spin superradiance. The pure spin superradiance is a purely self-organized 
process starting from an absolutely incoherent state when the average transverse mag- 
netization is strictly zero. The triggered spin superradiance is a process in which self- 
organization also plays an important role but whose beginning is triggered by an initial 
coherence imposed onto the spin system, that is by assuming that the mean transverse 
magnetization is not zero. 

The interpretation of pure spin superradiance cannot be based on the Bloch equa- 
tions because of the following. If the initial transverse magnetization is zero then, in 
the content of these equations, the relaxation of an inverted spin system can be due 
only to two reasons: either to spin-lattice interactions characterized by a relaxation 
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time Ti , or to thermal damping caused by the Nyquist noise of resonator. At very 
low temperature, typical of experiments [10-15] with polarized nuclear spins, the spin- 
lattice relaxation time 7\ is much longer than the dephasing time T 2 , therefore this 
mechanism cannot develop coherence. The resonator thermal damping, as is shown by 
Bloembergen and Pound [19], is negligibly small for macroscopic systems, the thermal 
relaxation time being proportional to the number of spins, N , in the sample, and 
so being much longer than not only T2 but even T\ . Thus, the resonator Nyquist 
noise can never produce the initial thermal relaxation. The radiation field in the coil 
does not provide a microscopic thermal relaxation mechanism, but the inhomogeneous 
internal, or local, fields are essential [19]. 

The Bloch equations cannot, in principle, describe the pure spin superradiance and, 
in general, any other relaxation regimes in which no initial coherence is imposed on the 
spin system. To treat all possible relaxation regimes for a nonequilibrium spin system, 
coupled with a resonator, it is necessary to take into account local spin fluctuations. 
This can be done by considering a microscopic model with realistic dipole interactions 
between nuclear spins. But, since the local fields are essential, we cannot invoke for 
a microscopic model a homogeneous approximation. The latter would immediately 
return us to the Bloch equations with the lost information on local spin fluctuations. 

If the number of spins, N , is not too large, say between 10 and 10 3 , then one 
can resort to a numerical solution of the corresponding evolution equations. Such a 
computer simulation, whose mathematical details can be found in [20], has been accom- 
plished [17] and confirmed the crucial importance of local spin fluctuations. These are 
sufficient for describing the pure spin superradiance, with no influence of the resonator 
thermal noise. 

Computer simulations, however, can give only a qualitative picture, as the number 
of spins involved is incomparably smaller than what one has in real samples with N 
of the order of 10 23 . In addition, such simulations provide no analytical formulae 
making it very difficult, if possible, to classify all possible relaxation regimes occurring 
when varying the numerous parameters of the system. 

The aim of the present paper is to untangle two mutually interrelated problems: 
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first, to formulate a method allowing an analytical solution for a system of nonlin- 
ear equations, with taking into account local fluctuating fields, as well as dynamic 
coupling and retardation effects; and second, to analyse various relaxation regimes of 
nonequilibrium nuclear magnets coupled with a resonator. 



The method to be presented here may be used not only for the particular problem 
discussed in the Introduction, but for a wide variety of evolution equations for different 
systems. In this section we will preserve the generality of the presentation. All neces- 
sary specifications related to the spin dynamics in nuclear magnets will be expounded 
in the following sections. To better understand the principal ideas of the method, it is 
convenient to divide it into several steps. 

1. Separation of variables 

Suppose that in the problem under consideration there is a set 



of small parameters. Depending on the way in which these parameters enter into the 
evolution equations, we may distinguish fast and slow variables. The terms describing 
local fluctuating fields can be treated as random, or stochastic, variables 



II. Method of Solution 



e = { £i \i = 1,2,...; 



Si\ < 1} 



(p = {ipi \i = 1,2,...; fi v } 



with a probability measure //^ . 



The fast variables 



u = {ui((p, t) \i — 1, 2, . . . ; t > 0} 



and slow variables 



s = { Sj (<p,t) |j = l,2,...; t>0} 



differ from each other by the properties of their evolution equations 



du 
~dt 



f(u,s,<f,t,e) 
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— = eg(u,s,<p,t,e 
whose right-hand sides are such that the limit 




(2) 



limf(u,s,ip,t,e) ^ 



(3) 



is not zero, while 



]imeg(u,s,<p,t,e) = 0. 



(4) 



Here and in what follows the matrix form of notation is used, according to which 
/ = {fi}> 9 = {di} j an d the product eg = c ij £ j9j} is to be understood as 
a column of linear combinations with coefficients c^- . All parameters, variables, 
functions, and coefficients can be complex except t > representing time. The limit 
e — > means that all £j — > . The right-hand sides of (1) and (2) can contain integral 
operators, provided that the limits (3) and (4) hold. For brevity, the dependence of 
the fast, u , and slow, s , variables on the parameters e is not explicitly written. 
Equations (1) and (2) are to be complimented by initial conditions 



The limiting properties (3) and (4) explain why the evolution equations of the form 
(1) correspond to fast variables, as compared to the evolution equations of the type (2) 
describing slow variables. 

The fast and slow variables are not necessarily simply defined for each given prob- 
lem, but the aim of this step is to introduce such variables by using the information on 
the existence of small parameters and by choosing the appropriate changes of variables, 
so that finally they could be distinguished in the above sense. 

2. Quasi-integrals of motion 

As far as the slow variables, by definition, vary with time much slower than the 
fast variables, the former may be considered as quasi-integrals of motion for the latter. 
Then we can try to solve the equations for fast variables under slow variables kept as 
fixed parameters. With the notation 



u(ip,0) = Uq, 



s(ip,0) = z Q . 



(5) 



u = X, 



s = 



(6) 
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where z is fixed, from (1) we have 

dX 

— = f(X,z, V ,t,e), (7) 

which defines 

X = X(z,<p,t). (8) 

The art of choosing variables is to get for (7) as simple equation as possible. In 
many cases this can be done so that (7), under fixed z , becomes a system of linear 
equations. The quasi-integrals of motion play here a role similar to the guiding centers 
in the guiding-center approach [21]. 

3. Method of averaging 

For the fast variable (8) we define the asymptotic period T by the condition 

Km\X(z,(p,t + T o )-X(z,(p,t)\=0. (9) 

If (9) gives several solutions for T , the smallest of them is to be taken. And if (9) 
has no solution for T , we put T — > oo . 

To find the time evolution of quasi-integrals of motion, we substitute (8) into the 
right-hand side of (2) and introduce the averaged function 



- r 1 f T o 

9(z,e)= — g(X(z,ip,t),z,ip,t,e) 
J |_io Jo 

Then the equation 



dt 



d/v (10) 



(ii) 

gives the sought time evolution. 

The foundation for this step is the Krylov-Bogolubov method of averaging [22,23]. 
The major difference in our case is that the Krylov-Bogolubov vector field (10) is de- 
fined as an average with respect to time and, in addition, with respect to the stochastic 
variable y? . 

4. Basic approximation 

The basic approximations for slow and fast variables are defined as follows. For the 
slow variables this is given by the solution 

z = z(t) (12) 



of equation (11), with the initial condition 



z(0) = z . (13) 

Substituting (12) into (8), we have 

x = x(<p,i)=X(z(t),<p,t) (14) 

for fast variables. The integration constant appearing when solving (7) is to be found 
from the initial condition 

x(ip,0)=u . (15) 

Note that (11) is,generally, a nonlinear equation, hence the basic approximations 
(12) and (14) take account of all nonlinearities essential for the considered dynamical 
process. 

5. Generalized expansion 

Corrections to the basic approximation can be found by using the generalized 
asymptotic expansion, 

oo 

u = x(<p, t) + J2 x n(v, t)e n , 

71=1 
OO 

s = z(t) + Y l z n (<p,t)e n , (16) 

71=1 

about (12) and (14). 

The right-hand sides of (1) and (2) are also to be expanded in a similar manner, as 



f(u, s, <p, t, e) = f(x, z, ip, t,e) + J2 fn{<P, t, e)e n . (17) 

71=1 

For example, in the first two orders we have 

h = Xlf'x + Zlfz, 

h = x 2 f x + Z2 f z + x lZl f^ z + \ (x\rL + z\r: z ) , 



where the notation 



d 

fx = faf( x > z > ^^,6); x = x(<p, t), z = z(t) 
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is used. 

The expansions (16) and (17) are to be substituted into the evolution equations (1) 
and (2). In doing this, we notice that, since because of (14) 

dx _ fdX\ fdX\ dz 
~dt ~ \~dt) z + \Jh) t Hi' 

then invoking (7) and (11), we get 

dx — 

— = f(x, z, <p, t, e) + e g (z, e)X' z (p, t), 

where 

X' z (ip,t) = ^X(z,ip,t)- z = z(t). 

Equating similar terms with respect to the power of e , we obtain the equations 
for the corrections of arbitrary order. It is important to stress that all these equations 
are linear, thus, there is no principal difficulty in solving them. To exemplify this, at 
the same time avoiding cumbersome formulae, let us think of e as of one parameter. 
Then for the first-order corrections we find the equations 

dx\ 

-^■ = f 1 (<p,t,e)-9(z,e)X' z { ( p,t), 

dz - 

—jj- = g(x,z,<p,t,e)- 9 (z,e). (18) 
The initial conditions, in compliance with (13) and (15), are 

Xl (<p,0) = 0, Zl (<p,0) = 0. (19) 

For all subsequent orders we have 

= fn(<P,t,e), 

^ = g n ( V ,t,e), (n>2), (20) 

with the initial conditions 

x n (^,0) = 0, z n (( P ,0)=0. (21) 
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The first of Eqs.(18) can be reduced to the form 



dt 



x 1 f x + A 1 -9X' z , 



in which 



Ax = h ~ xif x = zif' z . 
As we see, the equation for x\ is really linear, since 



Zl 



(<P,t) = J 



9(x,z,<p,t,e)- 9 (z,e) 



dt 



(22) 



immediately follows from the second of Eqs.(18). The solution for this linear equation 



is 



where 



Xl = e p e~ p [Ai-gX'Adt, 



(23) 



p = p(<p,t,e) = J f' x (<p,t,e)dt. 
For the second-order corrections, from (20), we find 

3 P | e~ p A 2 dt, z 2 = J gi dt 



x 2 = &' 



(24) 



with 



A 2 = h ~ *2f' x = z 2 f z + x lZl f^ z + - {x\f' x + z\f, 



Similarly, for the n -th order corrections we obtain the general formulae 



in which 



x n = e p J e- p A n dt, 
z n = j g n -idt, (n > 2), 



An — fn X n f x . 



(25) 



The simplicity of obtaining the higher-order corrections, satisfying linear equations, 
is a considerable advantage of the suggested generalized asymptotic expansion, as com- 
pared to the quiding-center approach [21] or averaging methods [22,23] in which each 
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subsequent approximation order invokes more and more complicated nonlinear equa- 
tions. Here we meet nonlinear equations only once, at the third step, when solving 
(11), which corresponds to the first-order averaging method. 

The use of the averaging method only in one step makes it possible, from one side, to 
include all essential nonlinearity into our basic approximation and, from another side, 
to define all corrections by simple formulae. The idea of dividing solutions onto their 
principal parts, including essential nonlinearities, and perturbative corrections, defined 
by linear equations, greatly helps in solving complicated nonlinear problems [24] . This 
idea, actually, goes back to the Struble technique [25,26] imployed for solving the 
Mathieu equation. Note that the nonlinear principal part could be also defined by other 
techniques known in the theory of singular perturbations [27], for instance, by using the 
methods of strained coordinates, multiple scales, nonlinear renormalizations, matched 
expansions, variation of parameters, and so on [28-30]. However, these methods, as is 
discueed in [31,32], are more ambiguous, more cumbersome, and less general than the 
method of averaging. 

Finally, we need to remember that, in our case, the solutions of nonlinear equations 
(1) and (2) contain the stochastic variable ip . As far as observable quantities should 
not depend on that variable, this means that the former are to be averaged with respect 
to the random ip with a given probability measure. The solutions themselves are not 
necessary such quantities that can be measured directly, but usually, the observables 
are some functions or functionals of these solutions. This especially concerns the fast 
variables, while the slow variables are often directly measurable. 



The system of nuclear spins can be modeled, as is accepted in the theory of nuclear 
magnetic resonance [2], by the Hamiltonian 



III. Nuclear Magnet 




(26) 



with the dipole interaction energy 




(27) 
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in which /i is a nuclear magneton, Si— {Sf, Sf, S?} is a spin operator, and 



r ij 



^ij — j T ij — T i 7" j i fij 

The total magnetic field 

B=H + H (28) 

consists of two parts, 

H = H e z , H= H 7 X , (29) 

the first is an external magnetic field H directed along the z -axis; the second, H , 
is a field of the coil of a resonance electric circuit, the coil axis being directed along the 
axis x . The sample is inserted into the coil. 
Introduce the interactions 

,2 



a ij = ~3 (l ~~ 3C0S 2 %) , 



3/i 2 

hj = —rj sin2 $ij exp(-i2ipij) , (30) 

Cij = --^-sin(2%) exp(-iipij), 

in which i?^- and ip^ are the spherical angles of . These interactions have the 
symmetry property 

(Xij djii bij ^jii *^ij Cji' (^"0 

Defining the ladder operators S~ and by the expressions 

$i = S'l ~ iSf i £>i~ = Si + iSf , 

s! = \{sr + s?), s?= l -(s--s+), (32) 

and using (30), we may cast the dipole interaction energy (27) into the form 
Hij = (yS^Sj — -S^Sj ^ + bijS^Sj~ + b i jS i Sj + 

+ 2 ( Cij St + dtjSr) S*. (33) 
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For the operators S { and S- , satisfying the commutation relations 



the Heisenberg equations of motion yield 

d N r / l \ 

fojisr = y: k,- (srs; + -s?ST )-2b ij sts;+ 

at jm l v z / 

+cy (SfS* - 2S?S?) + c*^ .S ; } - + ^5? (34) 



and 



j TV r 



+ ( Cij S+ - c^S-) S*} + (S; - St) - <R 7 i (5f - Ci) , (35) 

where (35) is supplemented by a term taking into account spin-lattice interactions 
leading to the longitudinal damping 7l , and Q being a stationary value of the spin 
z -component. The derivation of the spin-lattice term from microscopic spin-lattice 
interactions can be found in literature [1-3]. 

The initial state of the spin system is assumed to be nonequilibrium and character- 
ized by a statistical operator p(0) . So, the average spin 

(Si) = Trp(O) Si (t) = Tr/3(t) Si (0) 

is a function of time. The evolution equations for averages can be obtained by us- 
ing either the Liouville equation for the statistical operator p(t) or the Heisenberg 
equations of motion for operators. We prefer the latter way based on the Heisenberg 
equations (34) and (35). 

IV. Resonator Field 

The resonance electric circuit, coupled with the spin sample, is characterized by 
resistance R , inductance L and capacity C . The coil, in which the sample is 
immersed, has n turns of cross section A over a length / . The magnetic field 
inside the coil, 

* = Pa) 
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is formed by an electric current satisfying the Kirchhoff equation 

in which Ef is an electromotive force of external fields, if any, and of the thermal 
Nyquist noise; the magnetic flux 

An 

$ = —nA oV pM x (38) 
c 

is due to the x -component of the magnetization 

u N 

M * = -£E<<s?>; (39) 
iv i=i 

and the filling factor rj and spin density p are 

V N 

respectively. 

The resonance electric circuit will be called, for brevity, the resonator, and the 
internal coil field (36), the resonator field. For the latter, the Kirchhoff equation (37) 
can be rewritten as 



d H TT , dM x cE f 

_ + 2l3 H + J I H(r) iT = -4, w — + -±, (40) 



where 

1 / , Ann 2 An 



fhC \ cH 
is the resonator natural frequency, and 

R uj ( n — u ^ 
73 = 2L = 2Q V ~R" 

is the resonator damping. 

It is convenient to introduce the dimensionless resonator field 



h = f, (41) 



driving force 

cpEf 
nA hj 3 



f - (42) 
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and the dimensionless average magnetization 



M v 1 N 



— „ n%™> (43) 

in which v — x,y, z . Define the coupling constant 

«o = KV%—, (44) 

characterizing the strength of coupling between the spin system and resonator. Then 
the Kirchhoff equation (40) acquires the form 

^ + 2 l3 h + u? f h(r)dr = -4a ^ + isf. (45) 
at Jo at 

The resonator field h , as is seen from (45), can be induced by a driving force / and 
by moving, but not static, transverse magnetization. 

V. Average Magnetization 

The statistical averaging of a spin operator Sf = Sf(i) , with a = x, y, z , is given 

by 

(S?)=Trp(0)SHt)- (46) 

We shall use the notation 

u t = (ST), s t = (St). (47) 

The statistical operator p(0) in (46) defines the initial values of (47), that is, Mj(0) 
and Sj(0) . 

To obtain the evolution equations for the transverse, Ui , and longitudinal, Sj , 
magnetizations, we have to average the equations of motion (34) and (35), according 
to (46). The dipole interactions are of long-range type, therefore the double spin 
correlations can be decoupled in the mean-field approximation 

Although this decoupling is well justified for long-range forces [33], it has a deficiency 
that is important for nonequilibrium processes: it does not take into account the 
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attenuation due to spin-spin interactions. This attenuation appears in the higher- 
order corrections to the mean-field approximation. The derivation of the spin-spin 
damping 72 in the second-order perturbation theory can be found e.g. in ter Haar 
[1]. This damping has to be retained for a correct description of relaxation process, 
though 72 is much smaller than the Larmor frequency 

coo = ^ > 0. (48) 

At the same time the small second-order corrections to the oscillation frequency (48) 
can be neglected; alternatively, they can be included into the definition of ujq . The 
mean-field decoupling with corrections leading to the appearance of the spin-spin 
relaxation parameter 72 can be called the corrected mean-field approximation. Within 
the framework of this approximation, the averaging of (34) and (35) yields for the 
variables in (47) the equations 

at 

+ \ ( S i U 3 + ^ u i s j) - ZbijSiU* + Cij(UiU* - 2SiSj) + CijUiUjX (49) 



and 



i-^ = ^73/1 (ui - u*) + {dijiuiU* - u*Uj) + bijU*u* - b* j u i u j + 



dt 2 ,J y 11 h .... 

+(c i:j u* - c^u^sj} - (s - Q). (50) 
Introduce the arithmetic averages 

N -y N 

U = J7T, U ^ s =T?T, s i ( 51 ) 

i=l i=l 

for the transverse and longitudinal magnetizations, respectively, and also for a station- 
ary magnetization 

1 N 

Define 



N . . 
i=i 



1 N 

6 i = J^J2 ( ^ a H s i + C H U ] + 4 M i ) ( 53 ) 
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(Pi = E ( b H u j + c v s i) ( 54 ) 



which is a real quantity, and 

2 

which is complex. 

For the averages in (51), from (49) and (50), using the symmetry property (31), we 
find 

du 1 N 

i— = -(u + i>y 2 )u + 73/is + — E(<^ + <piSi) (55) 
d* iV i=1 

and 

ds 1 1 ^ 

^ = 2 7s/l ^ ~ ^ ~ ^ 7l ^ ~ ^ + 2iV S^i"* - ^ u <)- ( 56 ) 

The quantities (53) and (54) are local fluctuating fields [1], whose existence is 

due to the inhomogeneity of spin distribution. If one would resort to a homogeneous 

approximation, in which Uj and Sj do not depend on the index j , then 5i and tpi 

would be zero, since for the dipole interactions (30) we have 

N N N 

E a a - E K - E c v - 

i(#0 i(^) J(#<) 

when — > 00 and the spin sample is macroscopic in all three dimensions. The above 
sums can be nonzero if the number of spins is not high (N < 10) or if the sample has a 
specially prepared irregular shape. Then the nonzero values of these sums are defined by 
a nonuniformity in the space distribution of spins in the vicinity of the sample surface. 
Such a boundary nonuniformity for small, at least in one of dimensions, samples can 
lead to unisotropic effects in relaxation processes [16,34]. This kind of inhomogeneity 
of a sample inside a coil can be explicitly taken into account in the definition of the 
effective factor [19]. 

It is worth emphasizing that even when the spin sample is macroscopic and has a 
regular shape, so that the above sums over the dipole interactions (30) are nullified, 
nevertheless, the local fields (53) and (54) are nonzero if one does not invoke a uni- 
form approximation for the magnetizations Uj and Sj . The local nonuniformities 
contribute to the inhomogeneous dipole broadening [35]. What is the most important 
is that without taking into account such local fluctuating fields it is impossible, as 
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has been stressed by Bloembergen and Pound [19], to provide a correct description of 
relaxation in spin systems. 

At the same time, if (53) and (54) depend on the index i showing their local 
position, then the equations (55) and (56) are not closed, but for the case of N 
spins we need to deal with a system of 3iV equations defined in (49) and (50). For 
a macroscopic sample with N ~ 10 23 , to deal with such a number of nonlinear 
differential equations is a task that is not affordable even for a computer. 

A way out of this trouble is as follows. We may treat (53) and (54) as random 
fluctuating fields with a distribution given by a probability measure ^ . That is, we 
may put into correspondence to the local fields (53) and (54) stochastic fields 

{ip } <-> {5i}, {ip} <-> {ipi}, 

in which ip is real, representing the real 5i , and y? is complex representing the 
complex ipi . At the present stage an explicit form of the probability measure fi v is 
not important and will be considered later. 

With the stochastic representation of local fields in mind, equations (55) and (56) 
are reduced to 

du 

— = i(u - y? + ij2)u - i(^ 3 h + (fi)s (57) 

and 

ds % % 

-j t = g + <P) u* - - (ish + p*)u- 7l (s - C). (58) 

Since u is complex, the third equation, additional to (57) and (58), can be the equation 
for u* or for \u\ 2 . For the latter we have 

-^j\ u \ 2 = -27 2 |-u| 2 - i{^h + p)su* + 1(73/1 + (f*)su. (59) 

These equations are to be complimented by initial conditions 

u(0) = M , s(0) = zq. (60) 

Eqs.(57)-(59) for the magnetizations plus Eq.(45) for the resonator field form the 
basic system of equations permitting a correct description of relaxation processes for a 
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spin sample coupled with a resonator. The physical meaning of all terms in these equa- 
tions is quite transparent: The real random field tp shifts the oscillation frequency; 
and the term 73/1 + ip plays the role of an effective field acting on spins, h being the 
resonator field and ip , stochastic field caused by local fluctuations. If in (57)-(59) we 
would put (fo and f zero, then we would return to the Bloch equations; however 
the presence of these random fields, as is discussed above and will be demonstrated in 
what follows, provides a crucial relaxation mechanism. Note that the stochastic local 
fields interconnect the transverse and longitudinal components of magnetization, but 
do not change the absolute value of the latter whose time variation 

|(H 2 + S 2 ) = -27 2 H 2 -2 7 i S ( S -C) 

is caused only by the spin-spin dephasing collisions and spin-lattice interactions. 

If we would decide to invoke the adiabatic approximation, in the way one usually 
does, then we should put ^ — > in (57) which immediately results in the linear 
relation between h and u , that is, in the static approximation. However, as is 
discussed in the Introduction, such an approximation could be reasonable only at the 
final stage of relaxation, but cannot correctly describe transient phenomena. 

VI. Separation of Variables 

To solve the system of equations (57)-(59) and (45), we use the method developed 
in Sec. II. To this end, we need to separate fast from slow variables by defining the 
appropriate small parameters. Usually, the widths 71 and 72 are small as compared 
to u ; and 73 is small as compared to uj . The stochastic fields (fo and ip are 
also to be considered as small, since the corresponding local fields (53) and (54), as is 
evident from their definition, are of the order of the local dipole interactions, that is, 
of the order of 7* which is a part of the inhomogeneous dipole broadening; 7* being 
much smaller than u . Thus, there are four small parameters: 

(61) 

UJq LOq UJq LO 

An additional small parameter appears in the quasiresonance situation when the res- 
onator natural frequency is close to the Larmor frequency of spins. Then the detuning 
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from the resonance, A , gives another small parameter 

— <1 (A = u-uj ). (62) 

The quantities inverse to the corresponding widths define the characteristic times 

7\ = -, T 2 = -, T; = -, T 3 = -, (63) 
7i 72 7* 73 

among which Ti is the spin-lattice relaxation time; T 2 , spin-spin dephasing time; 

T 2 * , inhomogeneous dephasing time; T 3 , resonator ringing time. To be more cautions, 

it is worth noting that, in our case, the width 7* is due to local spin fluctuations which 

is only one of the possible mechanisms of inhomogeneous broadening. The latter arises 

also owing to crystalline defects, hyperfine interactions and other inhomogeneities [35] 

that are not included in our consideration. Therefore, here T 2 * is of the order of T 2 , 

both of them being related to dipole interactions, so 7* ~ 72 • The existence of the 

small parameters (61) means that the oscillation period 

2tt 

T = — < min{Ti, T 2 , T 2 *, T 3 } (64) 

is the shortest time as compared to the characteristic times (63). 

To check the properties (3) and (4), we have to take the limit in Eqs.(45) and (57)- 
(59) by putting zero all small parameters (61) and (62), and respectively, (po and if . 
This procedure yields the limits 



dt 

dh 



du 



ah, r 

— > — to 2 / h(r)dr — 2iaocu (u — u*), (65) 

dt Jo 



ds d . , 



which shows that u and h are to be treated as fast, while s and \u\ as slow 
variables. The first of limits in (65) also shows that the adiabatic approximation is not 
appropriate when u is not zero. 

At the next step we have to consider the slow variables as quasi-integrals of motion 
for fast variables. The corresponding equations (57) and (45), with the notation 

u = x — iy, s = z, (66) 
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where z is kept as a fixed parameter, can be written in the form 



dx 

— = - l2 x + u v y - <p 2 z, 



dy 



-u^x - 7 2 y + (73/1 + (pjz, 



(67) 




in which 



u} ip = u — ip 



(68) 



is the shifted frequency, and the stochastic field 



(69) 



is separated into its real and imaginary parts. The initial conditions to (67) are 



It is remarkable that the system of three integro-differential equations (67), under 
fixed z , is linear, thus can be solved exactly by imploying, e.g., the method of the 
Laplace transforms. Equivalently, differentiating the last of the equations in (67), we 
may convert (67) into a linear system of five ordinary differential equations, which is 
again exactly solvable by means of either the method of the Laplace transforms or the 
matrix methods. 

The exact solution of (67) is so cumbersome that it is not pleasure to write it 
down explicitly. Fortunately, we can simplify it by using the existence of the small 
parameters (61) and (62). Such a simplification can be done directly by, first, finding 
an exact solution of (67) and, second, performing some expansions in small parameters. 
However, this direct way is extremely tedious and does not provide an insight into the 
physics of the made simplifications. The same final result can be obtained in another 
way which is much less wearisome and more physically clear, and which is explained 
below. 

The formal solution of the last equation in (67) can be written as the sum 



x(0) = x , 



y(o) = yo, 



h(0) = 0. 



(70) 



h — h s + hf, 



(71) 
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in which the first term is a feedback field induced in the resonator by moving spins and 
the second term is a resonator field formed by driving forces. The resonator feedback 
field may be presented either as the convolution 

h s = -4a f ^x(t - r)W(r)dr (72) 
Jo at 

or as the Stieltjes integral 

h s = -4a [ Wit - r)dx(r), 

JO 

and the resonator forcing field is given by the convolution 

hf = ls fw{t-r)f{T)dr, (73) 
Jo 

where the transfer function is 

W{t) = (cos u 3 t - ^ sin u 3 t) (74) 

with 



^3 = 

The action of the resonator field (71) on the spin system involves, as follows from 
(67), the small parameter 73 . Neglecting this parameter reduces the first two equations 
in (67) to 

dx 

— = - l2 x + u v y - cp 2 z, 
dy 

— = -uj^x - >y 2 y + <piz. (75) 

The solution to (75) is 

x = (a cos uj^t + 60 sin uj^t) e~ 72< + —z, 

y = (b cos u^t — a sin u^i) e~ 72 * + —z, (76) 



0J<p 



where 



<Pi , ^2 
a = x z, b = y z. 
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Imploying (76) in (72) gives the feedback field 



h s — 

73 



dx ( tpi 

Oil— + Oi 2 U v X Z 



dt 



in which 



Oil 



OL 2 



«073(72 - 73) 



(72 - 7s) 2 + (A + y? ) 2 
"o73(A + <p ) 



72 



73; 



+ (A + ^o) 2 



3(72-73)* 



3(72-73)* 



(77) 



(78) 



If in the expression (77) we put oti = 0, a 2 = const , we return to the static- 
coupling approximation, while if we put oti = const, a 2 = , then we get the dynamic- 
coupling approximation [9]. However, in general, a.i = osi(t) and a 2 = ot 2 {t) are 
nonzero functions of time. The temporal dependence of the coupling functions in (78) 
portrays the retardation due to a gradual switching on of the coupling between the 
spins and resonator. Really, as is seen from (78), at the initial moment the coupling is 
absent 

ai (0) = 012(0) = 0. 

Using the first of the equations in (67) for (77) yields 
2 



h. 



73 



[(«i72 - Oi 2 oj tp )x - OLiOJ^y + (aiip 2 + at 2 (pi)< 



(79) 



Substituting (79) back into (67) reduces the system of three integro-differential equa- 
tions to the system of two ordinary differential equations 



dx 
~di 



-72X + uj^y - ip 2 z, 



dy 

— = -{My - 2a^ 2 z + 2a 2 u; lf z)x - (72 + 2a 1 u ip z)y+ 
+ (ifi + 2ai<p 2 z + 2a 2 (p 1 z + ^ 3 h f )z 

for the fast variables. 

VII. Fast Variables 

There is no problem in solving (80), which gives 

x = (ai cosfi^t + bi sinf^t) c~ Vipt + + Xf, 



(80) 
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y = (a 2 cos Qpt + b 2 sin Q^t) e~ T ^ + y v + yf, 
where the first parts describe the spin oscillations with the effective frequency 

1/2 



effective attenuation 



and coefficients 



(82) 



(83) 



a 1 =x -x v , a 2 = yo-y<p, 



^2 = [(1 + 2a 2 ^) + aizyo] + 



+ 



2ai z (^ + ai z^-)-(l + 2a 2 z)^- 

\ Win i Lib / ^ 'm 



y?2 



the terms 



* n 2 + r 2 



'1 + 2a 2 z)(fi - ( — - aiz J y3 2 



2/ v = 



+ 



1 + 2a 2 <2 + 2ot\Z 



. a;,, 



<P2 



(84) 



are originated by the local random fields; and the last terms 

x f =73 / Gi(t - T)h f {r)dT, 
Jo 

J// = 73/ G 2 (t - T)h f (r)dT 
Jo 

are due to the resonator forcing field; the Green functions being 



(85) 



G^t) = 2-^sin^t-e- r ^, 

^ Lip 

G 2 (t) = zcosrL v t ■ e~ Tipt - aizGi(t). 
In this way, the fast variable u , defined by Eq.(57), becomes 

U = U s + Up + Uf, 



(86) 



25 



where 

u s = {c x e iQjipt + c 2 e~ in ^ e~ r ^, 
u v = x v -iy v , Uf = x f -iy f , (87) 

and 

1 i 
ci = - (ai - 62) - - (61 + a 2 ) , 

1 i 
c 2 = - (ai + 6 2 ) + - (61 - a 2 ) ■ 

To find an explicit expression for Uf , induced by an electromotive force Ef , 
entering into the right-hand side of the Kirchhoff equation (37), we need to concretize 
the form of Ef . Accepting for the latter the standard expression 

Ef = E cos cut, (88) 

for the driving force (42) we have 



cfiE 
nA h^ 



f = f cosut; f Q = J±-J^. (89) 



Then the convolution (73), with the transfer function (74), gives 

h f = ^ (cos ut - ^ sin utj (l - e~ 73 ') . (90) 
Substituting the resonator forcing field (90) into (85), we get 

z/ = (/ie-* + /*e--*) (l-e-^), 

(f 2 e™ + fie-™) (l - e"^*) , (91) 



Vf 

where the coefficients are 



f, = foUJ ^ z (A +tT ) 



h = fi \i— - aiz) , 



and the effective detuning is 



k v = u-%. (92) 
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Therefore Uf in (87) becomes 

Uf = (d x e wt + d 2 e~ iu}t ) (l - e^ 3 ') 



(93) 



with the coefficients 



di = fi { 1 + — + ia x z 



do 



n ( i - ^ 



0J U 



+ ici\z 



Finally, the fast variable h , given by the sum (71), is composed of the terms (79) 
and (90) for which we have 



73 



i(a 1 + ia 2 )u* - i(a 1 - ia 2 )u H («i^2 + «2V 9 i) z 



UJu 



(94) 



and 



h f = ^(l - i^j (e^ + e~ i<JJt ) (l - e^ 3 ') . (95) 

The factors (1 - e" 73 *) in (90), (91), (93), and (95) describe the retardation in the 
interaction of the sample and resonator. 

VIII. Slow Variables 

At the next step of the method, displayed in Sec. II, we have to substitute the fast 
variables (86), (94), and (95) into the equations (58) and (59) for the slow variables 



s — z. 



\U\ = V, 



(96) 



averaging the right-hand sides of (58) and (59) over the asymptotic period of fast os- 
cillations and also over a distribution of stochastic fields characterized by a probability 
measure /i^ . The asymptotic period, according to the definition (9), is just (64). Let 
us denote the double averaging of a function F = F^t) , over the asymptotic period 
and over stochastic fields, as 

r T 



(97) 



Since ip is real and <p = ipi — iip 2 is complex, there are three independent real 
components of the stochastic fields, thence the differential measure dfj, v can be written 
as the product 

d^ = dii(ipo)dn{ipi)dn(ip 2 ). 
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It is customary to model the distribution of local dipole fields in spin systems by a 
Gaussian distribution [3,35]. Accepting this and assuming, for simplicity, that each 
distribution of ip u , with v — 0,1,2 , has the same width 7* , we get 

dfi((p v ) = — jj= exp \ -i ( — ' ' r 



Accomplishing the averaging (97), we will take into account the existence of the 
small parameters (61) and (62). The basic formulae that are met in the course of aver- 
aging the right-hand sides of (58) and (59) are assembled in the Appendix. Averaging 
the coupling functions in (78), we have 

_ // \\ f^A ^(72 -7s) 2 
a = {{ai)) = a 



Wq) (72 - 73) 2 + A 2 ' 

«"■» =o »(%) fa- ^A> - <98) 

The average effective frequency (82) and attenuation (83) are, respectively, 

n = ((n v )) = u (i + pz), 

r=((r ip ))= l2 + au z, (99) 

where an expansion in powers of the small parameters in (98) is used. 

To write the evolution equations for the slow variables (96) in a compact form, we 
shall use some notation. Introduce the effective coupling parameter 

_ ^0 vr(7 2 -7 3 ) 2 

g = a— = a — . 100 

72 \l2j (72-7a) 2 + A 

Define the damping 



ls = f -f* L + 2ny + -f^ M/?A - a l2 ) + y (aA + (3 l2 )] \ (101) 

oC^O [ ^ + 72 



appearing when calculating the correlator ((u s hf)) for the fields from (87) and (95), 
and also the attenuation 
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+ A^4~t| I (a " 27r/9)(A2 ~ ^ + 272 A(/9 + 27fa) l } (1 ° 2) 
resulting from the calculation of the correlator {{u/hf)) for the fields (93) and (95). 

Thus, the averaging of the right-hand sides of Eqs.(58) and (59), in compliance 
with (97), leads to the equations 

dz i /~\ 

— = g l2W - 7, - 7i (z - C) - IfZ, 

d ^ = -2 l2W -2{g l2 w- ls )z + 2 lf z 2 (103) 



for the slow variables, where 



w = v 2 -2e,z; e, = \ (104) 

The quantities (101) and (102) characterize the relaxation of the magnetization 
owing to the action of the resonator field (95) formed by driving force (89). Note that 
7 S = for the incoherent initial condition, when uq = xq — iyo — . The squared 
amplitude of the driving force (89), remembering (44), can be written as 

This shows that / ~ 1/y/N . Consequently, for the attenuations (101) and (102) 
we have 7 S ~ l/y/N and 7/ ~ 1/N . These values for a macroscopic sample with 
N r-j 10 23 should be negligibly small. 

In particular, if the electromotive force (88) corresponds to a resonance mode of 
the thermal Nyquist noise of the resonator, then [3] for its amplitude we have 

where ks is the Boltzmann constant and T , temperature. For uj in the radiofre- 
quency region, typical of spin systems, (106) simplifies to 

KM- (io7) 

Whence, for the amplitude in (105) we get 



2 



f^ = ~f~~^ {Nyquist noise). (108) 
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Substituting (108) into (101) and (102), we again come to the conclusion that these 
attenuations for a macroscopic sample are negligible. We shall exemplify this by nu- 
merical estimates in Sec.X. 

The conclusion that the radiation field of the coil does not provide a microscopic 
relaxation mechanism, so that 7 S and 7/ can be neglected in the equations for slow 
variables, is in complete agreement with the statement of Bloembergen and Pound [19] 
that a homogeneous magnetic field, such as exists in the coil, will never produce the 
initial thermal relaxation in a macroscopic sample. 

Let us acknowledge that 7 S and 7/ are negligibly small as compared to 72 . In 
addition, at low temperatures, characteristic of experiments [10-15], the spin-lattice 
damping is also much smaller than the spin-spin dephasing parameter. Thus, we have 

— <1, ^<1, — <1. (109) 

72 72 72 

Taking into consideration (109), the slow-variable equations in (103) can be contracted 
to 

dz 

Tt = m 

*jjL = -2 l2 w(\ + gz). (110) 

The equations in (110) can be solved exactly in the following way. Notice, that the 
effective attenuation (99), with notation (100), acquires the form 



T = l2 {l + gz). (Ill) 



Using (111) in (110), we obtain 



f = d i = - 2rro - (112) 



Differentiating the first equation in (112), we come to 

d 2 T dF 

+ 2T— = 0, 



dt 2 dt 

which yields 

dt 
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? + r 2 = 7o 2 , (us) 



where 70 is an integration constant. Eq.(113) is the Riccati equation whose solution 
is 

r = Totanh(^) (114) 

where t , having the meaning of a delay time, is another integration constant. From 
(111) and (114) we have 

z = ^tanhf^)-i (115) 
and from the first equation in (110) we find 

™=(^Vsech 2 (^y (116) 

The functions (115) and (116) are the exact solutions of (110). For the slow variable 
v , the relation (104) gives 

2 { 7o\ _„ U 2 ft — to 



v 2 — sech ^ + 2 £ ,z. (117) 

\9l2j V r / 

As is seen, r is an effective relaxation time. 

The integration constants 70 and t Q are to be found from the initial conditions 

z(0) = z , v{0)=v o . (118) 

From (115), (117) and (118) we obtain 

7 2 = r 2 + (( ? 7 2 ) 2 (t; 2 -2^o), 

r = 72(1 + gz ); 7ot = 1, (119) 

and the delay time 



+ r ° 1 
to = 2" ln 



70 - r. 



. (120) 

7o -r 1 

So, all constants in the solutions (115) and (117) for the slow variables are defined. 
The corresponding solutions for the fast variables are obtained by substituting (115) 
and (117) into the sums 

u = u s + + Uf, h = h s + hj, 
whose terms are given by (87), (84), (93), (94), and (95). 
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IX. Relaxation Regimes 

Depending on the initial conditions and system parameters, one can distinguish 
several qualitatively different relaxation regimes. The advantage of dealing with ana- 
lytical solutions, as compared to numerical solutions, is that there are explicit formulas 
allowing direct investigation. When the problem contains many parameters, as in the 
considered case, the detailed analysis of the solutions by varying the numerous param- 
eters becomes excessively laborous if not impossible. At the same time it may happen, 
that not all parameters are equally important, but only some of them or some their 
combinations. A striking example of this kind is presented by the problem considered 
here. Really, despite of great number of various parameters, characterizing the spin 
system coupled with a resonator, the solutions of evolution equations contain only 
several constants, the main of which is the effective coupling parameter (100). The 
general qualitative classification of different relaxation regimes can be done by varying 
only three quantities: the coupling parameter g , the initial polarization z , and the 
initial transverse magnetization v . The latter defines the level of initial coherence 
imposed on the system. 

First of all, one can easily observe that if there is neither initial polarization, nor 
initial coherence, than (110) has only the trivial solution 

z = v = (zq = v = 0). (121) 

Therefore, the necessary and sufficient condition for the existence of nontrivial solutions 
is a nonzero initial magnetization, 

ml = zl + vl> 0. (122) 

The relation between the effective relaxation time r and the spin-spin dephasing 
time T 2 depends on the value of gm . Namely, 

r ~T 2 (gm < 1), 

to < T 2 {gm > 1), (123) 
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which follows from (119) under the assumption that ge* <C 1 . The latter inequality 
is justified owing to the definition of in (104) as of a small parameter of second 
order with respect to (61). 

The delay time (120) can have either negative or positive sign depending on the 
value of gz Q : 

to < (gz > -1), 

to > (gz < -1). (124) 

If to < , then the maximum of the transverse magnetization (117) occurs at t = . 
In this case, since gz^ > — 1 , then To > , which means that the amplitude of the 
fast variable u decreases with time. When t > , then the maximum of (117), i.e. 
the maximum of coherence, occurs at t — t . In this situation, as far as gz < — 1 , we 
have T < , which leads, according to (112), to the increase of the amplitude of v . 
The negative sign of the attenuation T means that the system acts as a generator. 

Varying the quantities gz and gvo , we may distinguish seven qualitatively 
different relaxation regimes. 

1. Free induction: 

g\z \ < 1, < gv < 1; 

to < 0, to « T 2 . (125) 

This is the standard case of free nuclear induction, with the maximal coherence im- 
posed at t — and relaxation time T 2 . The coupling with a resonator plays no 
principal role. Note that the conditions of the upper line and lower line in (125) are not 
independent, but one line follows from another, in compliance with (123) and (124). 
However, we write down the relations between effective parameters, as well as those 
between characteristic times, to make the classification more physically transparent. 

2. Collective induction: 

gz > -1, gv > 1; 

t < 0, r < T 2 . (126) 

This case differs from the free induction by an essential role of the coupling with 
the resonator, which is sufficiently strong to develop collective effects leading to the 
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shortening of the relaxation time r . When gv ^> 1 ,then ro < T 2 . But, as in the 
previous case, the maximal coherence is that which is imposed at t — . 

3. Free relaxation: 

g\zo\ < 1, v = 0; 

to < 0, r « T 2 . (127) 

The initial polarization z and the coupling parameter g are not sufficiently high 
for the appearance of self-organized coherence. At the same time, there is no imposed 
coherence. The relaxation process is mainly incoherent being due to the local random 
fields. 

4. Collective relaxation: 

gz > 1, v = 0; 

to < 0, r < T 2 . (128) 

The difference with the previous case is that the positive initial polarization and the 
coupling parameter now are high, so that collective effects shorten the relaxation time. 
However, the initial state is close to a stationary one, and the change of v , being 
again due to the local fields, is too small to yield a noticeable coherence. 

5. Weak superradiance: 

-2 < gz < -1, v = 0; 

to > 0, r « T 2 . (129) 

The negative initial polarization corresponds to an inverted system. The value of this 
polarization and that of the coupling parameter g are sufficient to make the delay 
time positive and to develop a weak coherence, as a result of incipient self-organization. 
But the latter is not yet enough strong to shorten the relaxation time. 

6. Pure superradiance: 

gz < -2, v = 0; 
t > 0, r < T 2 . (130) 
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The system is prepared in a strongly nonequilibrium state with a high negative polar- 
ization. The coupling with a resonator is also strong. No initial coherence is imposed 
on the system. The coherence arises as a purely self-organized process started by local 
stochastic fields and developed owing to the resonator feedback field. 
7. Triggered superradiance: 

gz < -1, gv > 1; 

to > 0, r < T 2 . (131) 

The initial polarization is negative and the coupling with a resonator is strong enough, 
so that the collective behavior of spins, tight with each other through the feedback 
field, is important. But the relaxation is triggered by an imposed initial coherence. 
Therefore, this is a collective but not purely self-organized process. 

In this classification, three regimes, free induction, collective induction, and trig- 
gered superradiance are triggered by an initial coherence thrust upon spins, that is 
by setting v o ^ . Local random fields do not play an important role. Such kind 
of regimes can be described by the Bloch equations. Other four relaxation regimes, 
free relaxation, collective relaxation, weak superradiance, and pure superradiance, are 
initiated solely by local fields. No initial coherence is involved, i.e. v o = . The Bloch 
equations cannot treat these four regimes. 

Organizing the above classification, we separated qualitatively different relaxation 
types. As is clear, there can be intermediate kinds of relaxation in between these 
regimes. For example, the case when 

gz < -1, < gv < 1 

is between weak superradiance and triggered superradiance. In principle, everywhere 
in this classification the condition v — can be replaced by gv < 1 , to include 
the intermediate regimes. However, it seems reasonable to distinguish, first, different 
physical reasons causing different relaxation mechanisms. 

In the process of relaxation, the polarization (115), starting at z = zq , tends to 

z ~^(T 2 -t ) (i»t ). (132) 
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If the initial polarization zq is negative, then (132) shows that a noticeable polarization 
reversal to a positive value occurs for the case when r < T 2 , that is for pure and 
triggered superradiance; also, it may happen at collective induction, though then the 
initial polarization is not high. The highest initial polarization is needed for pure 
superradiance. The corresponding polarization threshold is twice as large as that for 
weak superradiance or triggered superradiance. Eq.(132) shows as well that there can 
be no essential reversal of polarization from positive to negative values. 

It is illustrative to consider more in detail two limiting situations, when the coupling 
of the spin system with a resonator is either weak or strong. Start with the weak 
coupling limit, g <C 1 . Then for the relaxation width and relaxation time, from (119), 
we get 

7o ~ 7 2 1 + gz + — (yl - 2e^ ) 



r ^T 2 
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For the delay time (120) we have 



^ln 



9 



The behavior of polarization is 



z - zo + |(«o ~~ 2e * z o) (l - e 



(133) 



(134) 



(135) 



When g\z \ < 1 and gv < 1 , we have the case of free induction (125), if v 7^ . 
And if v = , then we have free relaxation (127) with 

2 



9 f 



(l - e' 2 ^) . 



(136) 



z ~ z - ge*z 

The latter regime is entirely due to local fields, since if would be zero, then z ~ zo 
and there would be no relaxation. 

In the strong coupling limit, g ^> 1 , from (119) we find 



7o - 72 



ZQ 



ml 
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This, using the inequality e, C 1 , can be reduced to 

70 ~ gm 72 
T 2 



I _|_ z ° _ ( i _ 2:0 ' £ * 2;o 



To 



1 



gm\ 



+ 1 



3^o \ £*^o 



#™o7 m 



The delay time (120) takes the form 

ml(m - z )(gm - 1) - (gmp 5 - 2 )e*2;o 



*o - -g - ln 



moK + z )(gm + 1) - (gmg - z )£*2 
For the final polarization (132) at f >t we obtain 

1 



(137) 



(138) 



z ~ m — 



5 



1 h I 1 2 5 ,£ * 2; 



07715 



(139) 



These formulas for g^o > 1 , depending on the value of gzo , correspond either to 
collective induction (126) or to triggered superradiance (131). When v o = , we come, 
again depending on the value of gz , to collective relaxation, 

(128), weak superradiance (129) or pure superradiance (130). 

Note that if gzo < — 1 , then for any vq the maximal coherence is reached at 
t — to > , when 



9 



v(t ) ~ to . 



(140) 



To better emphasize the role of local fields, let us analyse the case when there is no 
initial coherence, that is 



TO = \z \, Vq = 0, 



(141) 



and g\z \ > 1 . Then 



7o - fi'Fo|72 



1 

1 + 

gzo 



T 2 




TO — . . 


fi-i + ( 


9 N 


gz \ 



l-J- * 

gzoj z 

L 3_\ £* 

gzo) z 



The delay time (120) becomes 



t 



n — — In 

2 



(\ z o\ ~ z ){g\z \ - 1) - (5-2:0 - l)e, 



(N + z )(g\z \ + 1) - (gz - l)e. 



(142) 



(143) 
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The final polarization (132) at t^$> t is 

1 



z ~ z - 



9\zo\ 



z \- z + (gz - 1)£*] . 



(144) 



Consider separately the cases of positive and negative initial polarizations. When 
the latter is positive, i.e. 

z = \z \, (145) 
then the delay time (143) and final polarization (144) are 



+ r o l 



(gz - !)£, 



2z (gz + 1) - (gz - l)e* 
z~z -(l — ) £* (i>i )- 



(146) 



Simplifying this for asymptotically large gz ^> 1 , and keeping in mind that <C 1 , 
we have 

T 2 



2^o 

2 ~ Z — 



r '- 
(t »t ). 



(147) 



Formulas (146) and (147) correspond to collective relaxation (128) due to local fields. 
Pass to the case of the negative initial polarization 



z o — —\ z o\ 



(148) 



Then, for the delay time (143) and final polarization (144) we find 



+ r ° i 
t ~—hi 



2N(#N - 1) + (g\zo\ + l)e* 



(g\zo\ + l)e* 



2 / 1 

Z ~ \Zq\ h 1H : r £„ 



(* »*o)- 



(149) 



This can describe weak superradiance (129) or pure superradiance (130). Under the 
inequalities g\z \ 1 and <C 1 the latter expressions change to 



+ r o l 
to In 



2^o 



T 2 

g\zo\ 



z ~ |2;o I h £* (£ < t ), 

g 



(150) 
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which corresponds to pure superradiance (130). The origin of this phenomenon is 
completely due to local fluctuating fields. 

An interesting question is: which part of dipole interaction is mainly responsible 
for starting the relaxation process in the regime of pure superradiance? Looking at 
Eqs.(57) and (58), we see that it is the random field (p which initiates the process, while 
(fo only shifts the oscillation frequency. The stochastic field <p represents the local 
fields (54), which are related to the terms bij and q,- of the dipole interactions (30). 
These terms are called nonsecular dipole interactions contrary to that is called 
the secular dipole interaction [2]. In this way, it is the nonsecular dipole interactions 
that originate an initial relaxation and, consequently, the pure spin superradiance. 

The obtained results make it possible to give one more justification for the term 
spin superradiance. For a system of N nuclei an effective number of radiators may 
be defined as 

iVe// - S ' 

where m is the initial magnetization introduced in (122) and S is nuclear spin. 
Averaging the power of current 

P v {t) = Rf = Np-h 2 , 
4a 

according to (97), we have 

try 2 

P(t) = ((P^t)))^N(a 2 + (3 2 )-^v 2 . 

a 

The average current power for a superradiant regime has a maximum at t — to > , 
where v(to) = m , is compliance with (140). Therefore, 

P(f ) ~ m 2 ~ N 2 eff . 

Also, as is seen from (137), the radiation time 

The situation when the radiation pulse is proportional to the number of radiators 
squared, and the radiation time is inversely proportional to this number, is character- 
istic of superradiance. 
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Note that the intensity of magnetodipole radiation I(t) , as a function of time, 
behaves similarly to the current power P(t) but contains a small factor making 
I(t) <C P(t) , so that P(t) is much easier to measure [17,19]. 

X. Numerical Estimates 

The aim of the present paper is not to discuss some particular experiments but 
rather to give the general picture of possible relaxation processes. Nevertheless, the 
general qualitative picture can be better understood if illustrated by quantitative esti- 
mates. For this purpose, let us accept the values of parameters typical of experiments 
[11-15] with proton-rich materials, such as propanadiol C 3 H 8 2 , butanol C^HgOH , 
and ammonia NH3 . Imploying the method of dynamic nuclear polarization, it is 
possible to polarize spins to a level of polarization reaching almost 100% . The sam- 
ples polarized in this way are good examples of metastable nuclear magnets. The 
lifetime of such metastable materials at low temperature is very long. This time, Ti , 
is related to spin-lattice relaxation time. The order of its magnitude is given by the 
relation 7\ ~ (a/A/) 2 T 2 , in which a ~ 10~ 8 cm is mean distance between spins, 
Al ~ 10~ 5 a ~ 10~ 13 cm is the coefficient of linear magnetostriction, and T 2 is the 
spin-spin relaxation time. Whence, Ti/T 2 ~ 10 10 . 

The spin-spin relaxation time is characterized by dipole interactions yielding T 2 ~ 
ha 3 //j 2 ~ 10~ 5 s . Consequently, 7\ ~ 10 5 s . The relaxation time T 2 * , related to local 
spin fluctuations, is also due to dipole interactions because of which T 2 * ~ 10~ 5 s . 

In principle, there exists another longitudinal relaxation time due to the interaction 
of spins through the common electromagnetic field formed under the magnetodipole 
spin radiation. This time, which will be denoted by T[ , to distinguish it from the 
spin-lattice relaxation time T\ , can be estimated as T{ ~ (A/a) 2 T 2 , where A 
is the radiation wavelength. For the external magnetic field H ~ 10 4 G , spins 
radiate in the radiofrequency region with u ~ 10 8 s~ 1 , thus with the wavelength 
A ~ 10 2 cm . This gives T[/T 2 ~ 10 20 or T[ ~ 10 15 s . As far as T[jT x ~ 10 10 , 
the longitudinal relaxation is practically due to the spin-lattice interactions only. The 
interaction through the radiation electromagnetic field is so weak, as compared to dipole 
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interactions, that it does not play any role. This drastically distinguishes spin systems 
from atomic and molecular ones exhibiting superradiance. In the latter systems, the 
effective interaction through the common radiation field is not only important but 
serves as the basic mechanism for the appearance of strong collective correlations and 
coherence. 

The resonator ringing time T 3 in the case of quasiresonance, when uj ~ uj ~ 
10 8 s~ 1 , and for the quality factor Q ~ 10 2 is T3 ~ 10 _6 s . The time of fast 
oscillations, defined in (64), is T ~ 10~ 8 s ; so it is really the shortest among other 
characteristic times. 

The damping parameters corresponding to the characteristic times in (63) are 71 ~ 
10 _5 s _1 , 72 ~ 10 5 s _1 , 7! ~ 10 5 s _1 , 73 ~ 10 6 s _1 . In this way, for the small parameters 
in (61) we have 

^~ht 13 , ^~io- 3 , ^i-io- 3 , ^~io- 2 . 

UJo UJo UJQ UJ 

The coupling constant (44), owing to the relations ^72 ~ /i 2 /a 3 and pa 3 = 1 , 
where p is the particle density, is ao ~ ^7772/73 ~ 10 _1 . The average coupling 
functions in (98) are a ~ 72/^0 ~ 10 -3 and j3 < 72/c^o ~ 10~ 3 . In the case of exact 
resonance, when A = , the latter is identically zero, (3 = 0. Thus, a and (3 are 
also small parameters. 

The maximal value of the effective coupling parameter (100) is of the order of n 2 . 
Therefore it varies in the interval < g oc 10 . 

Consider the dampings (101) and (102) caused by the action of the electromotive 
force corresponding to a resonance mode of the thermal Nyquist noise with the ampli- 
tude (106). The typical temperature in experiments [11-15] is T ~ 0.1K . As far as 
ksT ~ 10 _5 eV and Tiuj ~ 10~ 7 eV , we have huj/ksT ~ 10 -2 , hence the approxima- 
tion (107) is justified. Using ^73 ~ 10~ 9 eV , for the forcing-field amplitude (108) we 
find f ~ 10 2 /VN . Then, for the damping (101) we get 7, oc (10 5 /v / A 7 )s _1 . In the 
case of passive initial conditions, when x = yo = , the value of (101) is exactly zero, 
7 S = . Expression (102) yields 7/ ~ (10 7 /A^)s _1 . For a sample of about 1cm 3 
the number of protons is iV ~ 10 23 . Thence, the thermal-noise forcing field has the 
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amplitude / ~ 10~ 10 , so for the damping (101) and (102) we get ^ s < 10 _7 s _1 and 
7j ~ 10 _16 s _1 . These quantities are so much less than 72 that there is no any reason 
to keep them in the equations. This also concerns 71 . Really, the relations in (109) 
are 

^ ~ 10- 12 , ^ ~ 10- 21 , ^ ~ 10' 10 . 

72 72 72 

Therefore, the thermal Nyquist noise of a resonator has no influence on the spin dy- 
namics in a microscopic sample. 

One might ask a question: What should be the size of a sample on which the 
resonator thermal noise could produce a noticeable effect? This would happen if 7 S ~ 
72 , which gives iV ~ 1 , or when 7/ ~ 72 , from where N ~ 100 . For N > 100 the 
Nyquist noise is practically of no importance. 

The method of solving the equations, used in the present paper, makes it possible 
to take into account the retardation effects, related to the appearance of factors like 
(1 — e~ 73 *) . These effects are important for the correct description of relaxation 
processes. For example, the threshold of initial polarization for superradiance, weak or 
triggered, as follows from (129) and (131), is z ~ — 1/g . In percentage, for spin 1/2 
and g ~ 20 , this means that the superradiance threshold is —10% . Respectively, 
the threshold of pure superradiance, given in (130), is —20% . These values are 
in agreement with experiments [11-15]. While, if we would neglect the retardation 
replacing the factor (1 — e -73 *) by 1 , then for the superradiance threshold we would 
get — 72/tto^o = ~ 10~ 3 • In percentage, this makes —0.1% , which is 

unrealistically small. 

In the regime of pure spin superradiance, the characteristic times r and t can be 
estimated from (150). Since r ~ T 2 /g\z \ , taking g\z \ ~ 10 , we find the radiation 
time r ~ 10 _6 s . The local-field parameter, defined in (104), is ~ 10 -6 . Whence, 
for the delay time we obtain t ~ (3 -j- 5)r , that is t ~ 10~ 6 — 10~ 5 s . The reversed 
final polarization, according to (150), can reach 90% . Note that the problem of the fast 
polarization reversal of proton solid-state targets is of great practical importance for 
the study of scattering in high and intermediate energy physics [15]. The phenomenon 
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of spin superradiance can be used to achieve the desired fast repolarization. 
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Appendix 

Here we present the basic formulas for the averages defined in (97) and used in 
Sec. VIII when deriving the equations for slow variables. 

For the stochastic fields, with the Gaussian distribution in mind, we have 

«<A)» = «¥>i» = = 0, 

«¥>3» = «¥>?» = «¥>!» = ll 

Note that, instead of defining a particular distribution, we could postulate the above 
properties of random fields. 

In the following expressions the averaging (97) is accompanied by expansions in 
powers of small parameters (61): 

«e- r *'» * 1 - — , 

UJq 

«e( <A -»>*)> + — (iA- 72 ), 

C-**))) ~ -I — , 



{{e lult (l - 

((e iAt (l - e-^))) ~t3, 
(( e *(«+n*.-Hr*.)* (l - e ~ 73 '))) ~ -i— , 

2 47T 2 7| 



((e 2iw * (l -e" 73 *) 2 )) ~ (1 - 2ttz) 



, 7f 



2c 2 ' 

where and are given by (82) and (83), respectively, and Q with T are 
defined in (99). 

Emphasize the importance of the factor (1 — e -73 *) responsible for the retardation 
effects. 
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